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Abstract. We present an investigation of the structure of the emission line region in a sample of 12 single-peaked 
Active Galactic Nuclei (AGNs). Using the high resolution H/3 and Ha line profiles observed with the Isaac Newton 
Telescope (La Palma) we study the substructure in the lines (such as shoulders or bumps) which can indicate a 
disk or disk- like emission in Broad Line Regions (BLRs). Applying Gaussian analysis we found that both kinds 
of emission regions, BLR and NLR, are complex. In this sample the narrow [OIII] lines are composites of two 
components; NLRl which have random velocities from ~200 to 500 km/s and systematic velocities toward the blue 
from 20 to 350 km/s, and NLR2 with smaller random velocities (~100-200 km/s) and a redshift corresponding 
to the cosmological one. The BLR also have complex structure and we apply a two-component model assuming 
that the line wings originate in a very broad line region (VBLR) and the line core in an intermediate line region 
(ILR). The VBLR is assumed to be an accretion disk and the ILR a spherical emission region. The model fits 
very well the Ha and H/3 line profiles of the AGNs. 
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1. Introduction 

The concept of a disk geometry in the Broad Line Region 
(BLR) is very attractive because the most widely accepted 
model for Active Galactic Nuclei (AGNs) includes a super 
massive black hole fed by an accretion disk. The detec- 
tion and modeling of so me double-peaked Balmer lines 
has supported this idea llPerez et alJ 1198^ IChen et al " 



1989HChen fc HalDernlll989t lEracleous fc Haloernl 



Rodriguez- Ardila et alJ Tl996t IStorch i-Bergman n et alJ 
1997t iLivio fc Xii Il997t iHo et alJ 12000: Shields et alJ 



2000: S trateva et all 



200.lt 



Storchi-Bergmann et al.l 



Eracleous fc Haloernl 12003 1. However, the fraction of 
AGNs clearly showing double-peaked profiles is small 
and statistically insignificant. On the one hand, the 
existence of double-peaked lines should not be required as 
a necessary condition for the existence of a disk geometry 
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in BLRs. Even if the emission in a spectral line comes 
from a disk, the parameters of the disk (e.g. inclina- 
tion) can be such that one observes single-peaked lines 
llGhen fc Haloe rn 1989; Dumont fc Colhn-SouffrinI \l99(i 
iKollatschnvfc Bischoft.2002; .Kollatschnv.2003,i . Also, a 
Keplerian disk with disk wind can produce single-peaked 
broad emission lines as n ormally seen in most of the AGNs 
ijMurrav fc Ghianglll997iil . On the other hand, taking into 
account the complexity of emission line regions of AGNs 
(see e.g. Sulentic et al. 2000), one might expect that the 
broad emission lines are composed of radiation from two 
or more kinematically and physically different emission 
regions, i.e. that multiple BLR emission components 
with fundamentally different velocity distributions are 
present (see e.g. Romano et al. 1996). Consequently, 
one possibility could be that the emission of the disk 
is masked by the emis sion of another emission line 
region. Recently, IPopovid (|2003) investigated the physical 
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processes in BLRs using a Boltzmann-plot metliod, and 
found that probably 'physical conditions in regions which 
contribute to the line core and line wings are different'. 
This supports the idea that the broad optical lines 
originate in more than one emission region, i.e. that the 
Broad Line Emission Region is com plex and composed 
of at least two regions. Moreover, ICorbin fc BorosonI 
l)l996(l found that 'the difference between the Lya and 
H/3 full width at zero intensity (FWZI) values provides 
additional evidence of an optically thin very broad line 
region (VBLR) lying inside an intermediate line region 
(ILR) producing the profile cores'. Consequently we may 
expect that the VBLR can be formed in a disk or d i sk-like 
emission region. Moreover, recently IWang et al.l l)2003(l 
investigated the central engines of 37 radio-loud QSOs 
and found that their accretion rates suggest that most of 
the objects possess standard optically thick, geometrically 
thin accretion disks. In fact, recently ovic et al ]l2nniL 
120021 hnn^t) it was shown that broad emission lines of at 
least three AGNs (Akn 120, NGC 3516 and III Zw 2) can 
be well fitted with a model which has two components: 
(i) an accretion disk and (ii) a region with a geometry 
different from a disk. 

The aim of this paper is to test the validity of the 
two-component model of a BLR which contains an accre- 
tion disk and one additional emission region, i.e. to try 
to find evidence that suggests that the disk emission can 
contribute to the line emission even if the line profiles 
are single-peaked. To do this, we observed 12 AGNs in 
the Ha and H/3 wavelength region with the Isaac Newton 
Telescope (Sec. 2, throughout). High resolution spectra 
were analyzed first by Gaussian analysis (Sec. 3) and af- 
ter that we applied the two-component model for BLRs 
assuming that it is composed of a VBLR and an ILR (Sec. 
4). 

2. Observations and data reduction 

It is very important to notice that the two peaks produced 
by the disk may appear like two bumps in the blue and 
red parts of the H/3 and Ha line profiles. To find the sub- 
structure connected with disk emission one should obtain 
the spectral lines with a relatively high spectral resolu- 
tion and S/N ratio. We observed with the Isaac Newton 
Telescope (INT) 12 AGNs (see Table 1) which have been 
previously observed i n the X-ray band (Fe K g fine, see e.g. 
iNandra et J] l|l997l) : ISiilentic et aD l|l99al) ) and where, 
according to the X-ray results, one can expect that a disk 
geometry is present at least in the X-ray emitting region, 
i.e. that a disk exists whose signature might be observed 
in optical lines (emission of the outer part of the disk). 
The observed AGNs have no double-peaked Ha and H/3 
lines. 

The observations were performed with the 2.5 m INT 
at La Palma in the period of 21 - 25 of January 2002. 
The Intermediate Dispersion Spectrograph (IDS) and the 
235 camera (with chip EEVIO) in combination with the 
R1200Y (for the Ha wavelength region) and R1200B (for 
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Fig. 1. The estimate of the continuum (solid line) in the 
case of Mrk 1040 H/3 (top) and Ha (down). Dotted line 
the observed spectra and dashed line (below) the spectra 
after continuum subtraction. 

the H/3 wavelength region) gratings were used. The list of 
the observed AGNs with the coordinates, redshift, central 
wavelengths and the exposure times are given in Table 1. 
The seeing was around l".l and the slit width was 1". The 
spectral resolution was ^ 1 A. As one can see from Table 
1, we observed the Ha and H/3 wavelength line region for 
all galaxies, except Mrk 141 where only the Ha region was 
observed. Also, after calibration of the spectra, the H/3 line 
of Mrk 493 was too weak and the red wing of the 3C 273 
Ha line was too noisy, so for these two spectra we use the 
low resolution spectra observed with the HST (on Sep 4, 
1996 and Jan 31, 1999) with G400 and G750L gratings, 
respectively. 

CuNe and CuAr lamps were used for the wavelength 
calibration. Standard reduction procedures including fiat- 
fielding, wavelength calibration, spectral response, and sky 
subtraction were performed with the help of the IRAF 
software package. 

The software package DIPSO was used for reducing the 
level of the local continuum (by using the DIPSO routine 
'cdraw 1') fitted through the dots taken to be on the local 
continuum (See Fig. 1). 

The redshifts of the consi dered AGNs were taken from 
IVeron-Cettv fc VeronI tOO(t) . 

3. Line profile analysis 

To analyze the shape of the H/3 and Ha lines, we first fit- 
ted each line with the sum of Gaussian components. We 
used a minimalization routine to obtain the best fit pa- 
rameters. We also assumed that the narrow emission lines 
can be represented by one or more Gaussian components. 
In the fitting procedure, we looked for the minimal num- 
ber of Gaussian components needed to fit the lines. To 
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Object 


RA 


DEC 


z 


Central 


Date of 


Number of 


Exposures 


name 


h m s 


o ' " 




wavelength (A) 


observation 


spectra 


sec. 


Mrkl040 


2:28:14.3 


+31:18:40.4 


0.016652 


6742.3 


24-Jan-02 


4 


600 










4901.6 


25-Jan-02 


4 


1400 


3C120 


4:33:11.2 


+05:21:27.4 


0.033010 


6742.3 


24-Jan-02 


4 


450 










5027.9 


22-Jan-02 


3 


900 


NGC 3227 


10:23:29.6 


+19:52:15.2 


0.003839 


6742.3 


24-Jan-02 


4 


480 










4899 


25-Jan-02 


3 


600 


PG1116+215 


11:19:08.7 


+21:19:33.9 


0.175700 


5718.5 


23-Jan-02 


3 


1200 










7699.1 


24-Jan-02 


3 


500 


NGC4253 


12:18:26.1 


+29:48:57.5 


0.012929 


4475.6 


23-Jan-02 


4 


1200 










6742.3 


24-Jan-02 


8 


460, 500 


MrkllO 


9:25:11.3 


+52:17:29.4 


0.035291 


5027.9 


22-Jan-02 


3 


1200, 1500 










6742.3 


24-Jan-02 


4 


720 


Mrk 141 


10:19:10.0 


+63:58:36.6 


0.041673 


6744.8 


24-Jan-02 


4 


800 


3C273 


12:29:08.7 


+02:03:46.7 


0.158339 


5669.4 


25-Jan-02 


2 


550 










7699.1 


24-Jan-02 


3 


500 


Mrk817 


14:36:20.5 


+58:48:14.6 


0.031455 


4904.3 


25-Jan-02 


3 


550 










6742.3 


24-Jan-02 


3 


500 


Mrk493 


15:59:09.6 


+35:02:21.1 


0.031328 


6742.3 


24-Jan-02 


3 


360 


Mrk 841 


15:04:03.3 


+10:26:48.1 


0.036422 


6742.3 


24-Jan-02 


3 


450 










4901.6 


25-Jan-02 


6 


300 


PG1211+143 


12:14:19.4 


+14:03:44.0 


0.080900 


5027.9 


22-Jan-02 


3 


600 










6742.3 


24-Jan-02 


3 


500 



limit the number of free parameters in the fit we set some 
a priori constraints ((PoDOvic et alJl20ni[ 120021 1200.^) : 

- For the H/3 line in the fitting procedure we con- 
strained the Gaussian parameters as follows: 

(1) The three narrow Gaussians representing the 
[OIIIJAA4959, 5007 lines and the narrow H/3 component 
are fixed at the same redshift and the Gaussian widths 
are set proportional to their wavelengths. The Full Width 
at Half Maximum (FWHM) is connected with the width 
(W) of the Gaussian profile^ by FWHM = 2WVh[2; 

(2) We imposed the intens ity ratio of the two 
[OIII]AA4959,5007 lines as 1:3.03 llWiese et al.lll96fil) . 

(3) We included in the fit a red shelf Fe II template 
consisting of nine Fe II lin es belonging to the multiplets 
25, 36 and 42 jKoristaf 1992V We took the relative strength 
of these lines from lKoristal Ijl99^ and assumed that all Fe 
II lines originate in the same region, that is, that all of 
them have the same redshift, and widths proportional to 
their wavelengths. The template contains 9 Fe II lines from 
4855.5 A to 5018.4 A, and as one can see from Fig. 2 (e.g. 
Mrk 493 II/3), provides a sat isfactory fit to the iron lines 
in the H/3 wavele ngth region l|Koristall992l:IPOT ovic et al.l 
l200lLl2002ll2003l). In comna-rison with the list of Fe II lines 



given bv lSigut fc PradhanI l|2003l) . the Fe II template used 
here gives a better fit to the Fe II lines of the present sam- 
ple of A GNs. (4) We inclu ded the He I AA4921.93,5015.68 
A lines l|Veron et al.ll2002|) . assuming that they originate 
in the same region, i.e. that they have the same redshift 
and width. Here we should mention that the iron lines 



given as exp [-(AX/W)^] 



4923.9 A and 5018.4 A (included in the iron template) are 
very close to these two He I lines. Therefore it is difficult 
to find the contribution of these components. We can ob- 
tain a satisfactory fit of the line profiles without them, but 
in the case of Mrk 1040, Mrk 110, Mrk 493, Mrk 817 and 
NGC 4253 the He II lines can contribute to the red wing 
of H/3 (see Fig. 2, full lines in the red wing of H/3), espe- 
cially He IIA 4923.9 A, which tends to be stronger than 
He IIA 5018.4A. 

After a first attempt to fit the lines in the H/3 wave- 
length region, we recognized that the [QUI] lines were 
well fitted only in the central part, but the wings could 
not be appropriately fitted with assumption (2), i.e. the 
[OIIIJAA4959, 5007 lines show extended wings and cannot 
be properly fitted by only one Gaussian. We also notice 
that the wings are asymmetrical, with a shallower slope 
towards the blue. 

To find an acceptable fit of the [QUI] lines we first in- 
cluded two more Gaussians which followed constraint (2). 
However, we found that a satisfactory fit in some of the 
AGNs can be obtained only if the intensity ratio of the 
[OIII]A5007 and [OIII]A4959 is left as a free parameter. 
Since our main purpose is to find the best fit in order to 
subtract the satellite and narrow lines from the H/3 and 
Ha line profile, we therefore decided to leave as a free 
parameter the intensity ratio of the blue-shifted gaussian 
components of the [OIII]AA4959,50G7 lines. As a rule we 
find that the second Gaussian of the [OIII]AA4959, 5007 
lines is shifted to the blue and more broadened than the 
central narrow one (see Fig. 2). To check our results for 
the different components of the [QUI] lines, we subtracted 
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the contributions of Fe II and He II as well as the H/3 red 
wing from [OIII] using the DIPSO procedure for the con- 
tinuum (see Fig. 3). Then we fitted the [OIII] lines using 
two Gaussians and the ratio of the [OIII] lines 1:3.03. In 
this case we can fit the [OIII] lines of some of the AGNs in 
an adequate way (e.g. the [OIII] hues of 3C 120, Mrk 817, 
NGC 3227). There were small differences in the residue 
for Mrk 110, NGC 4253 and PG 1211. But for the re- 
maining two AGNs (Mrk 1040 and Mrk 841) we could 
obtain a satisfactory fit only with the intensity ratio of 
[OIII]A5007 and [OIII]A4959 left as a free parameter (see 
Fig. 4). Indeed, we have measured the fluxes of the [OIII] 
lines and found that in these two AGNs the intensity ratio 
tends to be smaller than 3.03. The results of the Gaussian 
fitting analysis for the [OIII]AA4959, 5007 lines are pre- 
sented in Fig. 4, and in Table 2. 

- For the Ha Une, we assumed that [Nil] AA6548, 6583 
and the Ha narrow component have the same redshift, 
and Gaussian widths proportional to their wavelengths. 
Taking into account that the two [NII]AA6548, 6583 lines 
belong to a transition within the same multiplet we as- 
sume an intensity ratio of 1:2.96 (see e.g. Wiese 1966). The 
fits of the Ha lines are presented in Fig. 5. One can ex- 
pect that other narrow lines in the optical spectra have the 
same shape as the [OIII] ones, but as the [Nil] AA6548, 6583 
lines are heavily blended we are not able to resolve their 
fine structure. Hence, each of them has been fitted by 
only one Gaussian. On the other h and, we include th e 
He IA6678.15 A and SiIIA6371 lines llVeronet al.ll2002^ . 
and each line has been represented by one Gaussian. Also, 
one can expect a contribution of Fe II lines in the Ha wave- 
length region, but on th e bases of the calculation given by 
ISigut fc Pradhai] l)2003|) there are no Fe II lines^ that can 
significantly contribute to the Ha line profile. Even if any 
lines are present their in tensity is probably negligible. It 
has be en shown (e.g. Ha loern fc O ke (19871): iLi'pari et aD 
(^2^) that AGNs with strong iron emission in the H/3 
wavelength region have a weak e mission of Fe I I in th e Ha 
band (see e.g. Figs. 1 and 2 in Hal pern fc Okd 1)19871) 1. By 
inspection of the NIST database'^ in the Ha region we 
found that a possible Fe II contribution may come from 
the permitted lines i(f - 3(f{^G)4s^G (Fe II lines 
6700.64 A and 6873.84 A ) and 3(f '^D2 - M^{^D)As'^D 
(Hues 6746.53 A , 6689.41 A and 6404.615 A ) as well as 
from FellA 6456.387. In the fitting procedure we included 
these lines assuming that they originate from the same re- 
gion, i.e. that they have the same width and redshift. As 
one can see from Fig. 5 (broad dashed lines), these lines, 
as well as the He I and Sill lines, make a negligible con- 
tribution to the Ha line profile, i.e. a satisfactory fit can 
be obtained without them. 



^ the nearest lines to Ha are Fe IIA6456.385 (near the blue 
wing) and Fe IIA7308.065 (near the red wing). 

^ http:/ /physics. nist. gov/PhysRefData/contents- 
atomic.html 



3.1. Discussion of the multi-Gaussian analysis 

In Fig. 6 we present the Gaussian widths of the different 
broad components versus their centroid velocities (relative 
to the systemic one) . The different components appear rel- 
atively well separated in this diagram, showing the consis- 
tency of the multi-Gaussian decomposition. By inspection 
of Figs. 2, 4 and 5 as well as of the diagrams in Figs. 6 and 
7, we can derive a number of conclusions concerning the 
broad line and the narrow line components, consequently 
concerning the BLR and the NLR. 

3.1.1. The Broad Emission Lines 

Considering only the broad line Gaussian components we 
can conclude: 

(i) that the H/3 and Ha line shapes of the considered 
AGNs are very complex, and usually cannot be described 
by one Gaussian, i.e. the Gaussian decomposition indi- 
cates a very complex kinematical structure of the BLRs. 

(ii) that the Gaussian decomposition indicates the ex- 
istence of a central broad component (part (II) of Fig. 6) 
with low random velocities (from 500 to 1500 km/s), and 
a redshift consistent with the systemic velocity (between 
± 500 km/s). 

(iii) that the red- and blue-shifted broad components 
(part (I) and (III) in Fig. 6) are often present in the fit. 
They tend to have higher random velocities with higher 
(positive or negative) shift. Such a near-symmetrical dis- 
tribution of the broad red- and blue-shifted components 
may indicate that in the considered AGN sample a frac- 
tion of the emission, at least the wings, originates in an 
accretion disk. 

(iv) The widths of the Fc II lines have different ran- 
dom velocities (Fig. 6) and a redshift consistent with the 
systemic velocity (±500 km/s). On the other hand, the Fe 
II random velocity (wpeii) tends to vary linearly with the 
central broad component random velocity (wc, see Fig. 7), 
with Wc < wpeii < 2wc except in the case of 3C 273 (the 
two points outside the region between these lines in Fig. 
7). This may indicate that a significant fraction of the Fe 
II emission is created in the region emitting the core of 
the Ha and H/3 lines. 

3.1.2. The Narrow Emission Lines 

The NLRs of the AGNs considered also show a complex 
structure, and we can clearly see at least two NLR regions: 

(i) the NLRl, which has an internal random velocity 
from ~ 200 to 500 km/s, and relative approaching veloc- 
ities from 20 to 350 km/s with respect to the systemic 
redshift of the observed galaxy; and 

(ii) the NLR2 which has an internal random velocity 
between 100 - 250 km/s, and a redshift equal to the sys- 
temic one of the corresponding object. 

On the other hand, the intensity ratio of the 
components from the NLR2 follows the line ratio 
(I(5007):I(4959)=3.03), while for the NLRl in 2 AGNs, 
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Table 2. The parameters of the Gaussian components of the [OIII] lines. Gaussian widths (W) and — z^Lf^ ~ 
znlri are given in km/s. The intensity ratio Rnlbi ~ ^5007/^^4959 of the 'blue component' obtained from the Gaussian 
fit, the measured flux ratio Rp = -^5007/^4959 and central intensity ratio of the 'blue' and 'central' Gaussian components 
(Ib/Ic) are also given (columns 2,3 and 6). The data for the objects denoted with stars should be taken with caution, 
because their [OIII] lines are too weak. 



Object 


Wnlri 


Rnlri 


Rf 


WnLR2 


Az 


Ib/Ic 


3cl20 


420 


3.0 


2.9±0.1 


130 


-110 


0.25 


3c273* 


950 


2.7 




497 


-200 


0.39 


Mrk 1040 


330 


1.8 


2.4±0.1 


190 


-340 


0.30 


Mrk 110 


300 


2.6 


3.0±0.1 


150 


-40 


0.17 


Mrk 817 


510 


3.0 


3.0±0.1 


145 


-300 


0.43 


Mrk 841 


215 


2.6 


2.8±0.1 


80 


-20 


0.45 


NGC 3227 


400 


3.0 


3.0±0.1 


140 


-120 


1.87 


NGC 4253 


340 


2.8 


3.0±0.1 


105 


-55 


0.55 


PG 1116* 


650 


1.8 




239 


-345 


0.88 


PG 1211 


270 


2.8 


3.0±0.1 


81 


-180 


1.8 



PG 1116 and Mrk 1040, this ratio is significantly smaller 
than 3.03 (see Table 2). Although a slight difference be- 
tween the observed and predicted intensity rati o of the 
[OIII] lines might exist ijStorev fc Zeippenl 12000 ^. this is 
certainly due to the bad S/N for PG1116, and in the case 
of Mrk 1040, to the high residues of the fit (see Figs. 2 
and 3). 

The difference in shifts and widths of 
[OIII]AA4959,5007 between these two NLRs indicate 
different kinematical and physical properties. Here we 
note that the parameters for the NLRl of 3C 273 and 
PG 1116+215 given in Table 2 should be taken with 
caution because of the low intensity of the [OIII] lines. 
We excluded these data in the discussion above. 

The clear tendency of the NLRl to have a blue-shifted 
systemic velocity (though in the case of Mrk 110 and Mrk 
841 the blue-shift is marginal and very close to the cosmo- 
logical redshift) supports the idea of a jet geometry of the 
NLR (see e.g. Dopita et al. 2003). In this case the receding 
jet component in the [OIII] lines might be obscured or ab- 
sorbed by the host galaxy, so one sees only the outflowing 
gas from the closer part of the jet. Many papers have been 
devoted to the radial velocity difference between narrow 
and broad lines (starting e.g. from Gaskell 1982), espe- 
cially the velocity difference between the [OIII] lines and 
H/3 that indicates a jet geometry of the NLR or of part 
of the NLR (see e.g. Bennert 2002, Zamanov 2002). On 
the other hand, 2D spectroscopy clearly shows that AGN 
NLRs have a complex structure (see e.g. Arribas 1997). 
The tendency for the [OIII] lines to be blue-shifted rela- 
tive to the up suggests that they are associated with a 
high-ionization outflow originating in these highly accret- 
ing sources (see e.g. Zamanov 2002). 

4. Two-component model for BLR 

Here we apply the two-component model in modeling the 
broad line shapes of 12 AGNs, where one component is 



the disk or disk-like reg ion. We start from the paper of 
ICorbin & BorosonI l)l996i'l . who investigated the combined 
ultraviolet and optical spectra of 48 QSOs and Seyfert 1 
galaxies in the redshift range 0.034-0.774. They found a 
statistically significant difference between the FWZI dis- 
tributions of the Lya and H/3 lines. The difference between 
the Lya and H/3 FWZI values provides additional evi- 
dence for an optically thin VBLR (which might be a disk 
or disk-like region) which contributes to the line wings. It 
is located inside an ILR which produces the profile cores. 
Also, they found relative weakness of the correlations be- 
tween the UV profile asymmetries and widths and those 
of H/3 line. This suggests a stratified structure of the BLR, 
consistent with the variability studies of Seyfert 1 galaxies 
(see e.g. Kollatschny 2003). The smaller average FWHM 
values of the UV lines compared to H/3 indicate that the 
ILR emission make a higher contribution to the UV lines, 
whereas in the Balmer lines the VBLR component is more 
dominant. This is also the case in well known AGNs with 
double-peaked Balmer lines, which usually show a single- 
peaked Lya line (see e.g. the case of Arp 102B, Halpern 
et al. 1996). 

The wings of the broad Ha emission line in the spectra 
of a large sa mple of AGNs (aroun d 100 spectra) were in- 
vestigated bv lRomano et alJ l)l996^ . They found a n indica- 
tion of multiple BLR emission. Moreover, recentlv lPopovid 
(2003) investigated the physical conditions in BLRs find- 
ing an indication that BLRs are complex and that physical 
conditions of the regions which contribute to the line core 
and the line wings are different. 

Though a two-component model can probably also be 
represented by other geometries, we choose the one with: 
a disk giving the wings of the lines, and a spherical com- 
ponent giving the core of the lines. 
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Fig. 2. Decomposition of the H/3 lines of the observed 
AGNs. The dashed hnes represent the observations and 
the solid lines show the profile obtained by Gaussian de- 
composition. The Gaussian components are presented at 
the bottom. The dashed complex lines at the bottom rep- 
resent the contribution of the Fe II, [OIII] and H/3 narrow 
lines 

4.1. Theoretical point 

For the disk we use the Keplerian relativistic model of 
Chen & Halpern (1989). The emissivity of the disk as a 
function of radius, R, is given by e = e^R^^ . Generally, 
when trying to fit the double-peaked line profiles by disk 
emission one should leave this index as a free parame- 
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Fig. 3. The estimated contribution of the H/3 wing and 
other satellite lines in the [OIII]AA4959,5007 line wave- 
length region (solid line) for Mrk 817. 
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Fig. 4. Decomposition of the [OIII] lines 

ter. However, we have to take into account two facts: (1) 
we have here single-peaked lines, i.e. the profile coming 
from the disk is not a priori well defined; (2) we will be 
using a two-component model which includes more pa- 
rameters than a disk-only model. We should therefore in- 
clude some constraints. Since the illumination is due to a 
point source radiating isotropically, located at the center 
of the disk, the flu x in the outer disk at diff erent radii 
should vary as ijEracleous k, HalDernlll994^ . We note 
here that this is indeed the way the incident flux varies, 
but n ot necessarily the way in which the lines respond 
to it llDumont fc Collin- Souffrii]lll990l iRokaki fc BoissonI 
ll999HEracleous fc Halpernll2003ll . However, the power in- 
dex p « 3 can be adopted as a reasonable prescription at 
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Fig. 5. The same as in Fig. 2. but for the Ha hne. The 
narrow dashed hnes at the bottom correspond to the nar- 
row Ha and [NH] Hnes, the broad dashed hnes correspond 
to the Fe H and He H contribution. 



least for Ha (|Eracleous fc Halr)ernll2003h . In our case, as 
one can see in Fig. 8, the shapes of Ha and H/3 in our 
sample are practically the same, so we can fit an averaged 
profile (see Sect. 4.2). Therefore we will start our fitting 
procedure by imposing the constraint p — but later we 
will change this parameter and use other constraints (see 
Sect. 4.2). 

We express the disk dimension in gravitational radii 
{Rg = GM/c^, G being the gravitational constant, M 
the mass of the central black hole, and c the velocity 
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Fig. 6. The widths (w) as a function of the inner shifts of 
the broad Gaussian components obtained for the Ha and 
H/3 lines of our sample of AGNs (open circles) and the Fe 
H template (full triangles) 
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Fig. 7. The Gaussian width of the Fe H template (wpeii) 
as a function of the central broad Gaussian component 
width iwc) of the Ha (open triangles) and H/3 (open cir- 
cles) line. The solid line corresponds to wpeii — wc and 
the dashed one to wpeii — 2 • wc- The points at the top, 
outside of the region between the two lines belong to 3C 
273. 

of light). The local broadening parameter (cr) and shift 
{zj) within the d i sk ha ve been taken into account as in 
IChen fc Halpernl l)l989tl . i.e. the 5 function has been re- 
placed by a Gaussian function: 



exp ■ 



(A-Ao-zO' 



2a2 



where zi and a are the local shift and the broadening 
parameter of the disk emission, respectively. 

On the other hand, we assume that the kinematics 
of the additional emission region can be described as the 
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Fig. 8. The comparison of H/3 and Ha (dashed Hnes) pro- 
files with an averaged one (sohd hne). 




Fig. 9. Two fits of 3C 273 with the two-component 
model. The disk parameters are: a) i = 14°, i?i„„=400 
=1420 Rg, Wd=1620 km/s, p=3.0 (Wg=1350 
km/s); b) i = 29°, i?„„=1250 Rg, Rout =15000 Rg, 
Wd=700 km/s, p=2.8 (Wg=1380 km/s) 



emission of a spherical region with an isotropic velocity 
distribution, i.e. with a local broadening wg and shift zq- 
Consequently, the emission line profile can be described 
by a Gaussian function. The whole line profile can be de- 
scribed by the relation: 

I{X) ^ IadW + IgW 

where IadW, IgW are the emissions of the relativistic 
accretion disk and of an additional region, respectively. 



4.2. Fitting procedure 

Before performing the fitting we 'cleaned' the spectra by 
subtracting (i) from H/3 the narrow H/3 line, the narrow 
[OIII] lines, the He I lines, and the Fe II template; (ii) 
from Ha the narrow Ha and [Nil] lines. Furthermore, we 
normalized the intensities of Ha and H/3 to unity and con- 
verted the wavelength into a velocity scale: A — > X = 
(A — Ao)/Ao. These conversions allowed us to compare the 
Ha and H/3 high resolution profile. We found that in the 
AGNs of the sample the Ha and H/3 have similar profiles 
(see Fig. 8). That concept supports the case that both hnes 
are formed in the same emission region. First we found an 
averaged line profile from the Ha and H/3 lines from each 
AGN (here we should note that for 3C 273 only the high 
resolution H/3, and for Mrk 493 and Mrk 141 only Ha were 
used). The averaged profile for each AGN was fitted with 
the above described two-component model. 

When a chi-square minimization including all the pa- 
rameters was attempted, it was found that the results are 
very dependent on the initial values given to the parame- 
ters. As mentioned above, the reason for this is that we 
apply a two-component model to single-peaked lines, so 
the number of free parameters is too large. To overcome 
this problem we have to use the additional constraint that 
the disc component fits the line wings, and the spherical 
component the line core. With this aim, we tried several 
values for the inclination. The starting point for the incli- 
nation can be the relative displacement of red- and blue- 
sh ifted Gaussia ns, Az. We use the empirical relation given 
bv lSmakI l|l98lh (see also Popovic 2002,2003). If we assume 
that a disk (or a disk-like) region exists, we can roughly 
estimate the parameters of the disk using the results of 
Gaussian analysis and the relationship (see Popovic 2002 
for more details) sini « Az\/2 ■ JSout, where i is the in- 
clination of the disk, and i?out is the outer radius given 
in gravitational radii. Taking into account that sini < 1, 
we can estimate the maximal outer radius. From Fig. 6 
one can estimate that 0.00033 (~1000 km/s) < Az < 
0.013 (~4000 km/s), thus the maximal outer radius may 
be in the interval ~ a few lO'^i?^ to ~ a few 10^i?g for the 
highest incli nation. On the other hand, from previous in- 
vestigations lichen et alJl989HEracleous fc HalperiJl994t 
lEracleousll2003l ) it was found that the outer radius of the 
disk has typical dimensions of ~ a few lO^i?^; in that case, 
the inclination for the sample is in the interval i ~ 5° — 25°. 

As we noted in Sect. 4.1, first we fixed p = 3, and chose 
different values of the inclination. Then we found the best 
fit by eye and after that we used a 'fixed' inclination 
such that other parameters can change more (with greater 
step) than the inclination'*. After that we performed a 
chi-square fitting of the parameters starting from suitable 
initial values. Note that this is not an 'orthodox' chi- 
square fitting procedure, but in this case the main point 
is not only to find the best fit, as we did with the Gaussian 

^ we use the SIMPLEX method for fitting that allows us 
to define the different magnitudes of 'refiections' (the steps to 
change the parameters). 
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fitting procedure, but also to try to explain the line profile 
with our two-component model. Consequently, we found 
that the fit of the wings of the BELs strongly restricts 
the value of the inner radius and the "local" broadening, 
i.e. the random velocity of emission gas in the disk. We 
found that the inclinations of the disk are relatively small, 
between 5° and 15°. One should have expected higher in- 
clinations in the sample (e.g. on average 30°). 

Therefore we decided to choose the inclination in the 
range i > 25°, leaving p as a free parameter, and look- 
ing for the maximal inclinations. We obtained reasonably 
good fits for the lines of more than half of the sample, for 
p < 3 and higher values of the inner and outer accretion 
disk radii (see Fig. 9ab). Note that l ow values of the emis- 
sivity index , p, were obtained also bv lEracleous Sz Haloernl 
l)l994|) and lEracleousI l)2003|) when fitting double peaked 
lines. Consequently, we did fitting tests for different fixed 
values of p in the range from 1.0 to 3.0. Also, in this case 
we can find acceptable fit for some AGNs. As an example, 
in Fig. 9 one can see two fits of the 3C 273 line with differ- 
ent parameters. The best fit corresponds to Fig. 9a, but 
the fit in Fig. 9b is also reasonable. From these fitting tests 
(without any constraints for the disk parameters) we are 
able only to give rough estimates of the disk parameters 
(see Table 3). We should mention here that for lines with 
the smallest asymmetry it was more difficult to estimate 
the parameters. Consequently, the estimated range of the 
parameters is higher (see Table 3, e.g. Mrk 110). 

5. Results and Discussion 

As one can see from Fig. 9, the line profiles can be well 
fitted with the two-component model, but some of the pa- 
rameters (e.g. the emissivity index, the inclination, the 
inner and outer radii) are not constrained. It is therefore 
not possible to find a unique solution for the model. For 
this, one should arbitrarily constrain at least one of the 
disk parameters^. However, it is obvious that in the AGNs 
of the sample, the shape of the line wings indicates radial 
motion, which may be caused by a disk-like geometry. In 
any case, the fitting tests described above allow us to ob- 
tain rough estimates of the kinematical parameters of the 
two-component model (Table 3). In Table 3 we present 
the estimated range of inclinations (i) and the minimal 
emissivity index (p"""), minimal and maximal value for 
shifts and widths of the Gaussian broadening term from 
the disk (2™'"'"^^''^ •y^mm.max-j^ shifts and widths of 
the Gaussian component (Z^Wg), and estimates for the 
minimal inner radius (R™J") and maximal outer radius 

l^out )■ 

Concerning the disk we can point out the following: (i) 
The maximal outer radius is in the range 10^ to 10^ Rg 
(ii) The local random velocities in the disk are different 
from object to object and they are in the broad range 
from ~ 300 km/s to 1700 km/s; so are the local redshifts 

^ e.g. the dimensions of the BLR obtained from reverberation 
studies might be used as a constraint. 
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Fig. 10. The random velocities of a spherical region (ILR) 
as a function of the local random disk velocities. The 
dashed fine represents the function Vjlr = Voisk km/s 



(zi ^ —780 to +450 km/s). (iii) The minimal inner radius 
of the emitting disk is in the range from ~ 100 to 600 
Kg, (iv) The inclinations are 5° < z < 50°, and in about 
half of the AGNs i < 30°. Such values of the inclinations 
support the idea that we more frequen tly observe the Sy 
1 type of AGNs at face-on inclinations ijRokaki fc BoissonI 
l1999il . 

Concerning the spherical emission region we can point 
out that: i) the redshifts are within 300 km/sec of the 
cosmological value; (ii) the random velocities in this region 
are also different for different objects, they are in the range 
from -400 to 1600 km/s. 

Particularly, it is interesting to see the correlation be- 
tween local broadening in the two regions presented in Fig. 
10. The local random velocities in the disk (averaged value 
from Table 3) are well correlated with those of the spher- 
ical region. This can be explained if this region originates 
from an accretion disk wind, which is created through sev- 
eral disturbances in the disk capable of producing shocks 
(e.g. Bondi-Hoyl flo w, stellar wind- wind collis ion, and tur- 
bulences) . Recently, iFromerth fc Mel ia* ('2001) described a 
scenario for the formation of a part of the BLR caused by 
shocks in the accretion disk; this may also create an ILR. 

We should mention that besides a disk (o r a disk-like 
region) or spiral shock way es within a disk ijChen et alJ 
Il989tlchen fc HalDernlll989l) . other geometries may cause 
the same kinds of substructure in the line profiles: i) emis- 
sion from the oppositely- directed sides of a bipolar outflow 
IjZheng et al .Il990tll99lh : ii) emission from a spherical sys- 
tem of clouds in randomly inclined Kepl erian orbits illumi- 
nated anisotropically from the center ijGoad fc Wanderl 
199(|); and iii) emiss ion from a binary black hole system 
(iGaskelll 1198,1 Il996 >). In any case one should consider a 
two-component model with an ILR contributing to the 
broad line cores and one additional emitting region con- 
tributing to the broad line wings. Recent investigations 
(see e.g. Wang et al. 2003, Eracleous & Halpern 2003) 
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Table 3. The parameters of the disk: zi is the shift and Wi = \/2a is the Gaussian broadening term from the disk 
which is a measure of the random velocity in the disk, Rinn is the inner radius, Rout is the outer radius, zq and Wg 
represent the parameters of the Gaussian component. 
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have shown that the disk geo metry for VBLR may be ac- 
cepted as a reality. Moreover . lEracleous fc HalperiJ l)2003|) 
found that the disk emission is more successful not only 
in explaining double-peaked line profiles but also in inter- 
preting the other spectroscopic properties of AGNs pre- 
senting these double-peaked Balmcr lines. 

6. Conclusions 

We observed 12 AGNs with INT in order to obtain high- 
resolution spectra of the H/3 and Ha lines that can be used 
for modeling BLRs. First we applied Gaussian analysis to 
the complex H/3 and Ha lines, from which we can con- 
clude: (i) both the BLR and the NLR are complex, (ii) 
the [OHI] narrow lines can be fitted satisfactory only with 
two Gaussians, one shifted toward blue, which may indi- 
cate the existence of an outflow, (iii) the different param- 
eters of these two Gaussians indicate that they originate 
in two NLR regions having different physical properties, 
(iv) the broad lines also show a complex structure, and 
they can be decomposed into three broad Gaussians - one 
red-, one blue-shifted and one at the systemic redshift. 

We adopt a two-component model, which comprises a 
VBLR and an ILR. We identify the VBLR with an accre- 
tion disk which contributes to the line wings. The cores of 
the lines originate in the ILR which is assumed to have a 
spherical geometry. This two-component model has been 
applied to the observed line profiles and we can conclude 
that : (i) The model can very well fit the observed line 
profiles, but it is very hard to obtain the disk parameters 
without imposing at least one constraint because of the 
large number of parameters and the lack of two peaks in 
the line profiles. They can be only roughly estimated us- 
ing fitting tests (see Table 3). (ii) The random velocities 
in the spherical emission region and the random veloci- 
ties in the disk are similar. This indicates that these two 
regions are linked through a common process, such as a 
wind produced by the disk. 

To find constraints for the model parameters further 
investigations are needed. 
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